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Abstract. 

We present an analysis of the dynamics of a nanomechanical resonator coupled to a 
superconducting single electron transistor (SSET) in the vicinity of the Josephson 
quasiparticle (JQP) and double Josephson quasiparticle (DJQP) resonances. For 
weak coupling and wide separation of dynamical timescales, we find that for either 
superconducting resonance the dynamics of the resonator is given by a Fokker-Planck 
equation, i.e., the SSET behaves effectively as an equilibrium heat bath, characterised 
by an effective temperature, which also damps the resonator and renormalizes its 
frequency. Depending on the gate and drain-source voltage bias points with respect 
to the superconducting resonance, the SSET can also give rise to an instability in 
the mechanical resonator marked by negative damping and temperature within the 
appropriate Fokker-Planck equation. Furthermore, sufficiently close to a resonance, 
we find that the Fokker-Planck description breaks down. We also point out that there 
is a close analogy between coupling a nanomechanical resonator to a SSET in the 
vicinity of the JQP resonance and Doppler cooling of atoms by means of lasers. 
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1. Introduction 

Nanomechanical single-electron transistors in which a mechanical resonator forms 
the voltage gate of the transistor constitute a new and interesting class of 
nanoelectromechanical system. The idea of coupling a nanomechanical resonator to 
the island of a single electron transistor (SET) as a mechanically compliant voltage gate 
was proposed as a way of measuring the displacement of a mechanical resonator with 
high precision [H El El El , since the conductance properties of the SET are extremely 
sensitive to the resonator motion. Indeed such devices have recently been used to 
measure the displacement of a nanomechanical resonator with almost quantum limited 
precision jB] . 

The sensitivity with which a SET can be used to measure the position of a 
nanomechanical resonator is ultimately limited by the back-action of the SET on the 
dynamics of the resonator. The back-action of electrons moving through a normal 
state SET gated by a nanomechanical resonator was studied recently [ZllHl El UHI and 
it was shown that, in the regime where the energy associated with the applied bias 
voltage is much larger than the resonator energy quanta, the SET electrons act on 
the nanomechanical resonator in a way which is closely analogous to an equilibrium 
thermal bath. In fact, the dynamics of the resonator can be described by a Fokker- 
Planck equation for a damped harmonic oscillator in contact with a thermal bath at 
a fixed temperature ^U]- Very similar results were obtained explicitly for a resonator 
coupled to a tunnel junction [HI [T2] and it was also shown by Clerk [in] that such 
behaviour is expected to be generic within the regime of linear response. 

In contrast to normal state SETs where the current arises solely from electron 
tunnelhng and cotunnelling processes [H], superconducting SETs (SSETs) can support 
a wide range of different electronic processes which contribute to the current including 
tunnelling or cotunnelling of quasiparticles, coherent tunnelling of Cooper pairs and 
even Andreev reflection [T^ ITHl [T71 ITHl IT^ 1^ . Furthermore, there exist a number of 
current resonances for particular values of the drain-source and gate voltages of the 
SSET where current is carried by a combination of different processes occuring at 
the source and drain junctions in turn. The best known (and most readily observed 
experimentally) current resonances for the SSET are the Josephson quasiparticle 
(JQP) and double Josephson quasiparticle (DJQP) cycles where transport occurs 
via a combination of coherent, resonant tunnelling of Cooper pairs and incoherent 
quasiparticle tunnelling US HnHmilHl 0120] . 

In this paper we analyse the back-action of a SSET on a nanomechanical resonator. 
In particular, we investigate the dynamics of a resonator coupled as a voltage-gate to a 
SSET which is tuned in the vicinity of the JQP or the DJQP resonance. We flnd that 
for both resonances the resonator can act as though it were coupled to a thermal bath 
with its dynamics described by a Fokker-Planck equation, as was found for the normal 
state SET. However, the magnitudes of the effective temperature and damping of the 
resonator in the vicinity of the JQP and DJQP resonances differ substantially, both from 
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each other and from those for the normal state case. A resonator coupled to a normal 
state SET has an effective temperature which is proportional to the drain-source voltage 
applied to the SET and is always damped. In contrast, the effective temperature of a 
resonator coupled to a SSET in the vicinity of the JQP or DJQP resonance is largely 
controlled by how far the SSET is biased from the resonance, rather than the magnitude 
of the drain-source voltage, and hence can easily be an order of magnitude smaller 
than the effective temperature for an analogous normal state SET. Furthermore, as the 
applied gate and drain-source voltages are adjusted to tune the SSET through a given 
JQP resonance, we find that the Fokker-Planck description breaks down sufficiently close 
to the resonance, while further from the resonance on the other side the Fokker-Planck 
description is restored once again, but with negative effective temperature and resonator 
damping constant implying the possibility of a dynamic instability. Very similar results 
for the SSET-resonator system have also been obtained independently using a different 
approach by Clerk and Bennett fIT\ . 

The dynamics of a nanomechanical resonator coupled to a SSET in the vicinity of 
the JQP resonance bears a striking resemblance to a number of other physical systems. 
In particular, the behaviour of the effective temperature of the resonator in the vicinity 
of the JQP resonance takes a very similar form to that of atoms undergoing Doppler 
cooling due to their interactions with laser light [221 1231 • Indeed, the minimum effective 
temperature of both a resonator in the vicinity of the JQP resonance and Doppler-cooled 
atoms are given by a decay rate: the quasiparticle decay rate for the SSET-resonator 
system and the decay rate of the excited state for the atoms. 

This paper is organized as follows. In section 2, we introduce a master equation 
describing the coupled statistical dynamics of the SSET-resonator system in the vicinity 
of the JQP resonance (a similar master equation for the SSET-resonator system in the 
vicinity of the DJQP resonance is described in the appendix). We then show that the 
master equation can be well-approximated by a Fokker-Planck equation under conditions 
of weak coupling and wide separation of SSET and oscillator dynamics timescales. In 
section 3 we present analytic and numerical calculations of the SSET-induced damping, 
frequency renormalization and effective temperature in the vicinity of the JQP and 
DJQP resonances. In section 4, we discuss our results and the analogy between 
the SSET-resonator device and other physical systems, before we finally present our 
conclusions. 

2. Master Equation description for the JQP resonance 

In this section we obtain a master equation for the SSET-resonator system in the vicinity 
of the JQP resonance and show that the dynamics of the resonator can be described 
by a Fokker-Planck equation. The same approach can also be used to derive analogous 
results for the DJQP resonance, details of which are given in the appendix. 

The model circuit that we consider is shown in figure ^ The SSET consists of 
a small superconducting island and two superconducting leads weakly-linked to the 
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Figure 1. Model circuit of the SSET-resonator system. 



island via tunnel junctions with capacitances Cj; a drain-source bias voltage Vds is 
applied to the leads. The nanomechanical resonator is treated as a single-mode harmonic 
oscillator with frequency uj and mass m. The metallized resonator is located adjacent the 
SSET island, forming a gate with capacitance Cg{x), which depends on the resonator's 
displacement x; a gate voltage Vg is applied to the resonator. In the experiments of 
Refs. El; the coupling between the SSET island and the resonator is typically very 
weak so that the displacement of the resonator from its equilibrium position is much 
less than the separation d between that equilibrium position and the SSET island itself. 
Hence if we also assume a parallel-plate geometry for simplicity, the gate capacitance 



can be approximated by jH] : Cg 



X] 



Cg{l—x/d), implying linear coupling between the 



SSET and the resonator [T^ . 

The central island of the SSET is taken to be sufficiently small that its charging 
energy = e^/2(2Cj + Cg) ~ A ^ k^T, where A is the superconducting gap and T is 
the temperature of the quasiparticles in the leads. Hence, the number of charge states 
accessible to the island is severely restricted. The Josephson coupling between the leads 
and the island is Ej = hA/{8e'^Rj) -C Ec, where Rj the resistance of the junctions |19j . 
Depending on the exact value of the polarization charge induced on the SSET island by 
the resonator gate, Ng = {CgVg + CjV^s)/e, certain quasiparticle and resonant Cooper 
pair tunnelling processes can become energetically favourable leading to a number of 
possible current carrying regimes such as the JQP and DJQP cycles. At sufficiently 
large drain-source voltages, and for relatively low junction resistances, it is also possible 
for current to flow via higher order processes such as quasiparticle co-tunnelling, but we 
will neglect such effects in what follows. 



Dynamics of a nanomechanical resonator 

(a) JQP Cycle 



5 



N+2 




o 

N+l 



(b) DJQP Cycle 




JN+2 
N+1 



o 



HI' 



Figure 2. Schematic illustration of the JQP and DJQP cycles, (a) For the JQP 
cycle, Josephson tunnelling involving a Cooper pair occurs between drain and island 
electrodes, increasing the island electron number N by 2, followed by two, subsequent 
quasiparticle tunnel decay processes into the source electrode, decreasing the electron 
island number by 2 [^. (b) For the DJQP cycle, Josephson tunnelling of a Cooper pair 
occurs at both junctions, with a Cooper pair tunnelling event through a given junction 
alternating in turn with a quasiparticle tunnelling event at the opposite junction jl9| . 



The details of the specific electronic processes which occur at the JQP and DJQP 
resonances are illustrated schematically in figure |21 Close to the JQP resonance, 
Cooper pairs tunnel between the right (with reference to the circuit in Figure Q), drain 
electrode and island, while electron quasiparticles tunnel out from the island to the 
left, source electrode. Alternatively, Cooper pair tunnelling can occur between the left 
source electrode and island, while electron quasiparticles tunnel in from the right drain 
electrode to the island. Which of these two JQP cycles takes place depends on the gate 
and drain-source voltage biases. We shall consider biases such that only the former cycle 
occurs (i.e., that illustrated in figure |2K)- 

We seek a master equation that describes the dynamics of the island charge state of 
the SSET and the position-velocity state of the resonator's center-of-mass valid in the 
vicinity of the JQP resonance. Master equations for the island charges in normal state 
and superconducting SETs have been derived by a number of groups using essentially the 
same technique jTHIIHI- Starting with the full (time-dependent) Schrodinger equation 
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for the system, one traces over the microscopic electronic degrees of freedom making 
use of a second-order Born approximation (that treats the quasiparticle tunneUing rates 
between the island and the leads as a small expansion parameter), followed by a long- 
time limit Markov approximation (that treats the response of the electrons in the leads to 
a tunnelling event as being instantaneous), to arrive at a master equation for a reduced 
density matrix represented in the basis of total number of electrons, N, on the island. We 
generalise this approach to include the resonator. Again starting from a fully quantum 
Hamiltonian, we follow the same procedure of the Born and Markov approximations 
and trace over the microscopic degrees of freedom. However, we also assume that 
the resonator does not evolve at all on the time-scale of the quasiparticle tunnelling 
processes. Essentially this means that we are treating the resonator as a classical 
oscillator|. Finally we take the Wigner transform j23] of the resulting equations to 
obtain the desired master equation which can be thought of as providing a semiclassical 
description of the coupled dynamics j^. 

For the JQP process (see figure Et), the island electron number can be A^, + 1 
or A^ + 2, with the A^ and A^ + 2 number states linked by coherent Cooper pair 
tunneUing. Hence the associated set of coupled semiclassical master equations has 
diagonal components Pn{x, v, t), pn+i{x, v, t), pAr+2(x, v, t), and off-diagonal component 
Pn,n+2{x, V, t) = p*j^^2 Af(^) where x and v are the position and velocity coordinates, 
respectively, of the oscillator. In our semiclassical description, the sum pN{x,v,t) + 
Pn+i{x, V, t) +Pn+2{x, V, t) is the probability density Piio{x, v, t) of finding the oscillator 
at the point in phase space (x, v) at time t, while the integral J dxdvpN{x, v, t) gives the 
probability P/v(t) that the island electron number is A^ at time t, with the probability 
conservation P/v + -P/v+i + Pn+2 = 1- Explicitly, the semiclassical master equations take 
the form 



I In deriving the quasiparticle tunnelling terms we effectively treat the resonator as a classical oscillator 
by making an adiabatic approximation, i.e. we assume that it does not move on the time-scale of the 
tunnelling processes, an approximation which was also used in Ref. 0. The other terms in the master 
equations which arise from the coherent evolution of the resonator and SSET charge are not affected 
by the adiabatic approximation. As we shall see later on, the position dependence of the quasiparticle 
transition rates in fact do not play an important role in affecting the resonator dynamics near the JQP 
resonance so we expect our master equations to provide a description of the resonator dynamics close to 
the JQP resonance which is essentially the same as that which would be obtained from a fully quantum 
derivation. 



Pn+2 



Pn 



Pn+1 




iPN+2,N — pN,N+2) 
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Pn) 



(1) 



where Xg = 2EcCgVg/ (emu'^d), which is the distance between equihbrium positions of 
the oscillator with and + 1 electrons on the island, parametrises the strength of 
the electro-mechanical coupling. The quasiparticle tunnel rates are given by [2^1 



is the normalized quasiparticle density of states, with B(-) the stepfunction, and 
n(e) = 1/ [1 + exp(e/A;B7')]- In the above master equations, the typical position 
coordinate is assumed sufficiently small that the tunnel rates can be expanded to first 
order in the position coordinate, with T'{E) = dV/dE. The quantities E]^^2,n+i and 
En+i,n are the energies gained by an electron when it tunnels from island to source 
electrode, the island number changing from + 2 to 1 and 1 to A^, respectively. 
The quantity Ej^+2,n is the energy gained by a Cooper pair when it tunnels from island 
to drain electrode. These energies are as follows: 



The JQP resonance condition is Enj^2,n = which is satisfied for A^^ = A^ + 1. 
Furthermore, the bias voltage must be sufficiently large for the quasiparticle processes 
to be allowed, enabhng the JQP cycle. For superconductors at zero temperature, this 
translates into the requirement that eV^s > 2A + E^. Notice also that in the absence of 
the mechanical oscillator, we recover the standard master equations for the SSET about 
the JQP resonance [T^ . 

It is convenient to express the master equation in terms of dimensionless 
coordinates, since in dimensionless form the essential parameters governing the dynamics 
are more clearly expressed. Rewriting the time coordinate in units of the tunnelling time, 
^tunnel = ^Rj/Vds, the positiou coordinate in units of Xg, and the velocity coordinate in 
units of Xs/rtunnei, the master equations take the form 





where 




(3) 



En+2,n+i = - 2E,{Ng -N- 3/2) + eV^s 
En+i,n = - 2E,{Ng -N- 1/2) + eV^s 
En+2,n = -^E,iNg-N-l). 



(4) 



PN+2 



Pn 
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dpN+i dpN+i 

— V- 



dv dx 



PN+2 



Pn,N+2 



PN+1 

+ 2nir (^En+2,n + Pn,n+2 

~ o r(-EAr+2,Af+l) + f '(_E'Ar+2,Ar+i)KX 



+ ivrej {pn+2 - Pn) 



Pn,n+2- 



(5) 



The definitions of tlie various dimensionless parameters are as follows: eno = ^^Ttunnoi 
is the ratio of the SET quasiparticle tunnelling time to the oscillator period, ej = 
Ttunnei/TRabi = {^Rj /Vds) {E j / k) = A/(8eVrfs) is the ratio of the quasiparticle tunnelling 
time to the Cooper pair Rabi oscillation period, r = Rj /{h/e^) is the ratio of the tunnel 
junction resistance to the quantum of resistance, and k = muj'^x1/{eVds) characterizes 
the coupling strength between the oscillator and the SET. The dimensionless tunnelling 
rate is 



/oo 
d~tp{l)p{i + E)n{l)[l-n{l + E)l 
-oo 



where now 



P{e) 



AO ^ 



e2 - A2 



(6) 



(7) 



and n{e) = l/[l+exp(eeKi,/A;BT)], with e = e/{eVds), E = E/{eVds), and A = A/{eVds). 
Note that for typical nanomechanical-SSETs [3 El, we have eno ^ej^l, /t^l, and 
r > 1. Also, the dimensionless quasiparticle tunnel rates F, and their gradients F', are 
generally of order unity. 

Our goal is to obtain a description of the dynamics of the resonator, decoupled 
from the details of the electronic degrees of freedom. One very direct way to obtain the 
resonator dynamics is to solve numerically the above master equation for the oscillator 
probability density Piio{x,v,t). Another direction is to take advantage of the typical 
conditions of weak coupling (k -C 1) and wide separation of timescales (eno ^ 1) to 
derive from the above master equation a much simpler, approximate effective equation 
for the oscillator probability density Piio{x, v, t) (i.e. a reduced master equation for the 
oscillator alone) which can then be easily solved. We have used both approaches. In the 
remainder of this section we describe how the reduced master equation is obtained and 
show that it is nothing other than the Fokker-Planck equation. Later in section 3 we 
compare the results of this approach with direct numerical integrations of the original 
set of master equations. 

We begin our derivation of the reduced master equation of the resonator by rewriting 
the full set of the master equations (0) in the following, concise 5x5 matrix operator 
form: 



v = {no + v) V, 



(8) 
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where 



V 



( pN+2ix,V,t) \ 

Im pN,N+2{x,V,t) 
\ Re PN,N+2ix,V,t) J 



Ti-o 



"RO 



d d 

X- V- 



( 



+ 



v 



dv dx 

rA?"+2,Ar+i 


^N+2,N+1 

nej 




(9) 







-Trej 






riv+i,Ar 
— rjv+i,Ar 





-27rej 
27rej 


rAf+2,jV+l 

2 

-2'KrEN+2 



N 







2nrEN+2,N 

rjv+2,]v+i 



with X denoting the 5x5 identity matrix and we use the shorthand notation 
r(£^jv+2,Ar+i) = rjv+2,iv+i (we have dropped the tilde for convenience, understanding 
that all quantities are in dimensionless form). The operator describing the interaction 
between the SET and oscillator is V = Vi + V2 



KxUi + e\Q^U2 where 



/ -r 



N+2,N+1 





r', 



^+2,^+1 


















r' 



iV+l,7V 
^ N+l,N 










■ jV+2,Af+l 
2 

— 47rr 






47rr 



r' 



N+2,N+\ 



(10) 



and 



Uo 



/ 1 + AP 



















-1 + AP 



















AP 



















AP 







V 











AP 


/ 



(11) 



where we have used the shorthand notation AP = (P/v) — (PAr+2). Note that we have 
redefined the position coordinate such that its origin coincides with the steady-state 
value (x) = — (A^ + 1) + (Pn) — (-P/v+2), where the steady state island occupation 
probabilities are taken to be those for the limit k, 0, 



{Pi 



N+2/ — 



(r^+2,iv+i/2)2 + (27rr£;^+2,iv)^ + {Trej^ (2 + ^^^) 



(Pn) 



N+2,N+1 



N+2) 



'^N+1,N 

1-{Pn+i)-{Pn+2), (12) 

an approach which is valid for sufficiently weak coupling. The advantage of working 
with this redefined position coordinate will become apparent shortly. 
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Equation (jSI) resembles the time-dependent Schrodinger equation. The 
'Hamiltonian operator' Tig gives the free, decoupled evolution of the independent 
oscillator and SSET systems, while the operator V = Vi + V2 describes the interaction 
between the two systems with Vi giving the dependence of the Cooper-pair and 
quasiparticle tunnelling rates on the oscillator position and V2 giving the SSET island 
number dependence of the electrostatic force acting on the oscillator. 

Given the close resemblance of equation (jH)) to the Schrodinger equation, we can 
apply approximation techniques developed for open quantum systems, in particular the 
self-consistent Born approximation (SCBA), followed by the Markov approximation. 
Applying the SCBA as described in section 3.1 of Ref. [27j, assuming weak coupling 
between the oscillator and SSET, k <^ 1, we obtain the following approximate expression 
for the reduced master equation probability distribution Pho(3;, f , t) of the oscillator E 

Pho(x, V, t) = HuoPnoix, v, t) + e^«°*TrsET [V(t)PsET(0)] e-'^^°'Puo{x, v, t) 



+ 



dt' e«Ho<TrgET [V(t)PsET(0)] TrgET [V(t')^SET(0)] e-^Ho*PHo(x, v, t) 
dt' e«H°*TrsET [V(t)V(t')^SET(0)]e-^«°*PHo(a;,t;,t), 



where T^ho 
and V{t) = 



av 



— V 



dx 



^Hot 



(13) 



is the Hamiltonian operator for the free harmonic oscillator 
is in the interaction picture. The initial, t = probability 

PHo(a;,t',0)PsET(0), where 



distribution is taken to be a product state: V{0) 

( PiV+2(0) \ 

Pset(O) = Pjv+i(0) 

Im Pn,n+2{^) 
\ Re P7v,7v+2(0) ) 

Note that the above SCBA step which gives the oscillator master equation p3|) and also 
the Markov approximation applied below should not be confused with the Born-Markov 
approximation described earlier in section 2 which gives the starting oscillator-SSET 
master equation these two approximation steps rely on distinct weak coupling and 
timescale conditions. 

The SCBA approach was applied to obtain a Fokker-Planck equation for a resonator 
coupled to a normal state SET in Ref. jTO] and our derivation for the SSET case follows 
the same route. We can use a Markov approximation and evaluate the integrals in ()13p 
for t ^ 00 since typically cho ^ 1 and we are only interested in the oscillator dynamics 
on timescales of order the mechanical period and longer, t > e^. Furthermore, using 
the redefined position coordinate, we find that the second and third terms on the right 
hand side of equation (fT!^ drop out and we eventually obtain 



(9P 



HO 



dt 



odd 2 



d 



21 



X 



Cfv) 



d_ 
' dv 



d 



d 



+^HO o . l-l a . + 1-2 o_ 



dv 



dx 



HO, 



(14) 
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where the C coefficients are defined as follows: 



and 



Ci^ = lim [ drTr[W,(t)W,(t-r)PsET(0)] (15) 

Ci' = lim [ dTrTT[Ui{t)mt~T)VsET{0)]. (16) 

In the Markovian limit which we have used, these coefficients do not depend on the 
initial state 'Pset(O) of the SET, just as we would expect. Furthermore, the d'^P/dvdx 
term in equation (fT^ (called the 'anomalous diffusion' term in Ref. |2Z|) is of order e 
smaller than the diffusion term when time is expressed in units of the oscillator period; 
it should have only a small effect on timescales of order the mechanical period or longer. 
Re-expressing equation (|Tl|l in terms of dimensionful coordinates and dropping the 
anomalous diffusion term, we obtain 

dPno 



dt 



u^x— -V— + ^ V I ^set^bTset ^ 
^ dv dx dv m dv"^ 



HO; 



(17) 



where the renormalized oscillator frequency is 



UR = ^^l + t^Cf (18) 
the damping rate is 

7SET = -neu;Cl\ (19) 
and the effective SET temperature is 

fceTsET = -eVds^. (20) 

'-'2 

Equation (fT7|) . which has the form of a particular class of Fokker-Planck equation 
known as the Klein-Kramers equation j2H|; describes the Brownian motion of a harmonic 
oscillator interacting with a thermal bath. The oscillator experiences a net damping 
force, due to the interaction with the SSET, [the third term on the right-hand-side of 
equation (fTTj) ] and an accompanying Gaussian distributed thermal fluctuating force [the 
fourth term on the right-hand-side of ()17p . called the 'diffusion' term]. 

Note that our derivation of the Fokker-Planck equation (fT7|) for the oscillator does 
not in fact depend on the speciflcs of the SET interacting with it. As long as the original 
master equation has the form V = {Ho + V)V with V = V1+V2 = KxUi+e^Q^U2 where 
operators lAi{2) involve only the SET parameters, then the same effective thermal bath 
description ()14j) of the SET results. The speciflcs of the SET enter in the dimensionless 
C coefficients deflned in equations ()15|1 and Thus, for the example of a normal 

state SET with island electron number fluctuating between the values N and iV + 1, the 
coefficients take the values = = —1 and CP = {Pm){Pn+i) [El CHI- Furthermore, 
the same approach can be applied to extract the relevant resonator dynamics for the 
DJQP cycle (see section 3.3 and the appendix). 
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3. Results 

The coefficients Ci and C2 , given by equations (fTH|l and (fTH|l . together with the Fokker- 
Planck equation (fT^ . describe the dynamics of a resonator coupled to a SSET. In 
this section we evaluate equations ()15|1 and (|T6|l numerically to obtain the effective 
temperature, renormalized frequency and damping of the resonator close to the JQP 
and DJQP resonances. In order to get a better picture of the underlying physics, and 
also to understand the limitations of the SCBA approach, we compare these numerical 
results for the JQP approach with analytical approximations and numerical results 
obtained by direct integration of the master equations [given by equation ©]• 

We begin this section by presenting approximate analytical results derived for the 
JQP resonance as they provide a useful framework within which to understand the 
resonator dynamics. We then present our numerical results for the JQP and DJQP 
resonances. 



3.1. Analytical approximations for the JQP resonance 

It is possible to derive analytical approximations to the C coefficients as series expansions 
in ej. Neglecting also the F' terms in the interaction operator Ui, since they are an order 
of magnitude smaller than the Airr terms in Ui, we find to leading non- vanishing order 
in ej: 

^21 _ (47r^rej)^£'7v+2,Ar {^n+2,n+i + '^^n+i,n) ^2]^^ 

^N+l,N [(rAr+2,7V+l/2)2 + (27rrEAr+2,iv)2]^ ' 

^21 ^ (47r^rej)^£'jv+2,AfrAf+2,Af+i 

['^'^N+l,N + (rAr+2,Af+l/2)^ + rAr+l,ArrAr+2,Ar+l + (27Tr_Ejv+2,Ar)^] , 



^22 ^ (7rej)^rAr+2,Af+l ^23) 

' r^^,^^ [(r^+2,JV+i/2)2 + (27rrE^+2,^)2]2 

^ [2r^_^]^^^ + (r7v+2,Af+l/2)^ + rAr+l,ArrAr+2,Ar+l + (2vrr£'7v+2,Ar)^] • 

The approximations to the renormalized frequency, damping rate and effective SET 
temperature then follow by substituting (j2H), (j221), and into equations (fTH|l . (jTHl), 
and ()2U|). In particular, for the SET temperature, we have: 

ksTs^ = -eV,.^ = J^^;^^^ . (24) 

where rAr+2,Ar4.i and E^^2,n are in their original, dimensional form [see equations ^ 
and (jl}]. Note that Tset does not depend on ej to leading order, while 7set and the 
frequency renormalization are 0(ej). 

Although these results are simple and intuitive, they are in principle only valid 
when vrej ^ 1 a condition which is by no means always satisfied for SSETs in practice. 
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3.2. Numerical results for the JQP resonance 

Two distinct sets of numerical calculations are carried out. We carry out integrations 
of the original master equations and within the framework of the Born-Markov 
approximation (0), we evaluate the expressions for the coefficients [equations (fT3j) and 
(fT?)|l ] numerically. 

For the numerical integration of the master equations jH], we obtain the full 
evolution of the resonator probability distribution (with initial state chosen to be 
Gaussian) from which we determine the evolution of the average position (with respect to 
the fixed point value), {x(t)). We then obtain values for uji and 7set by fitting (x(t)) to 
the equation of motion of a damped harmonic oscillator with a renormalized frequency. 
Finally, we use equations (fTHj) and (fT^ to infer values of the corresponding coefficients. 
In obtaining these results we concentrate on the typical experimentally accessible regime 
K, e -C 1. As a consequence, the numerical integrations take a prohibitively long time 
to reach a steady-state and so we do not extract values of Tset from the integrations. 

Figures El and ID show our numerical calculations of the coefficients Cf^, 
and TgET, with the corresponding analytical approximations (PT|) . and 
shown for comparison. Notice that while the numerical results for Cf^ and were 
obtained from numerical integrations of both the master equations (without Born- 
Markov approximation) and from equations (fTH|l and (fTH|l (with the Born-Markov 
approximation), the numerical results for Tset were only obtained using the latter 
technique. In the numerics, values for the system parameters were chosen that are typical 
of those found in current devices 0121 • In particular, the value of ej = A/ (8eV^s) = 1/16 
is kept relatively small. This value corresponds to choosing eVds = 2A, the correct order 
of magnitude to enable the JQP and DJQP cycles. Nevertheless, it is clear from figure El 
that our analytical approximations differ substantially from the numerics close to the 
centre of the resonance. In contrast, figure HI shows that there is excellent agreement 
between the analytical and numerical calculations of Tset- It is not clear why the 
agreement is so good, but the most obvious explanation is that the cancellation of the 
ej dependent terms that occurs in our approximate expression for the ratio CP/C|^ must 
extend beyond second order in ej. 

Figure El shows schematically the Vds and Vg bias ranges in relation to the JQP 
resonance lines for which the curves in figures El and |31 are obtained. Note that, because 
the quasiparticle tunnel rates and their gradients have been approximated as constants 
equal to unity for simplicity, the JQP curves in figures El and IH do not depend on the 
Vds bias choice. With the dependences of the tunnel rates on the energies Ep^+2,N+i 
and En+i,n properly taken into account [see equations and (jHl)], one finds that the 
maximum and minimum values of the Cf^ and curves increase in magnitude as Vds 
decreases towards the onset for the JQP cycle at eVds = 2A + Ec- In contrast, as we 
shall see below the DJQP curves depend strongly on the Vds bias choice even when the 
quasiparticle tunnel rates are approximated as constants. 

The results obtained from the numerical integrations of the master equations agree 
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Bom-Markov 
Analytical Approx. 
1^0.01 




E /eV 

N+2.N ds 



Figure 3. Comparison of numerical calculations of the coefficients and 
with the analytical approximations [equations H21(l and ()22|l]. The full curves come 
from numerical evaluations of equations H15|l and (|16|l . whereas the points come from 
numerical integrations of the master equations ©• For the numerics we choose 
ej = 1/16, r = 1 and set the quasiparticle tunnel rates and their gradients to unity. 
In addition, when integrating the master equations we set e = 0.1 and vary the value 
of K as shown. 



well with those obtained within the Born-Markov approximation. As expected, the 
agreement between the two improves as the magnitude of k is reduced, in accord with 
our use of the approximation k -C 1 in deriving the reduced master equation. 

In addition, the numerical integration of the master equations gives us one more 
piece of information that we could not have obtained from our calculations within the 
Born-Markov approximation. Very close to the resonance we find that the evolution of 
{x{t)) no longer matches that of a damped harmonic oscillator. §. Thus, it seems that 

§ The points plotted in figure |2| are all obtained from a fit to the behaviour of a damped harmonic 
oscillator. Close to the resonance where this fit is not readily obtained we were not able to extract 
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Figure 4. Comparison of numerical calculations of Tset with the analytical 
approximation, equation (|24|l . The parameters are the same as those used for figure|31 




Figure 5. Schematic map showing the location of a selection of neighbouring JQP 
resonances (solid black lines) and DJQP resonances (solid black circles). The bias 
ranges for the C^^, and Tset plots are indicated by the sohd red hues. 
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E /eV 

N+2,N ds 



Figure 6. Numerical calculations of the coefficients Cf^ and in the vicinity of the 
DJQP resonance. The different values of Ec/eVds correspond to points at different 
distances below, (a) and (c), and above, (b) and (d), the centre of the DJQP resonance 
(see figurel^J. As before, we choose ej — 1/16, r — \ and set the quasiparticle tunnel 
rates and their gradients to unity. 



very close to the resonance the action of the SSET on the resonator is not analogous 
to a thermal bath. Although we have not carried out a systematic investigation of the 
region in which deviations from thermal bath-like behaviour occur, we do find that the 
width of the region (in terms of Ej^+2,n) broadens with increasing k. This suggests that 
the weak coupling (Born) approximation we use to derive the Fokker-Planck equation 
may break down for mujQxl > En+2,n (i-e., k > En+2,n)- 

3.3. DJQP resonance 

The coefficients and calculated within the Born-Markov approximation are 
shown in figure El while the corresponding variation in Tset is shown in figure [7| To 
facilitate comparisons with the JQP resonance we have chosen E^j^2,n and eVds/Ec 
as independent variables, in terms of which the other energy detuning is Ej^j^i^^-i = 

values for 7set and hence no corresponding data points are plotted for C21 
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Figure 7. Numerical calculation of Tset (a) below and (b) above the D JQP resonance. 
The parameters are the same as those used for figure 



En+2,n — 4:{Ec/eVds) + 2. The relation between the range of SET bias points covered in 
figures ini and [7| are illustrated schematically in figure El 

The DJQP resonance is much more complex than the JQP one and our results by 
no means cover the whole range of the relevant parameter space. However, by comparing 
figures El and ini it is clear that the magnitude of C|^, and hence 'Jset, is generally much 
larger for the DJQP resonance than for the JQP one. As we shall discuss below, this 
implies that in an experiment the effects of the back-action due to the SSET on the 
resonator dynamics will be stronger close to the DJQP resonance. 

By analogy with the JQP resonance, it seems inevitable that there will be regions 
very close to the DJQP resonances where our Born-Markov approximation will fail to 
capture fully the physics of the system. However, because of the additional complexity 
of the DJQP cycle, we have not as yet performed the necessary numerical integrations 
of the master equations required to determine exactly where and how the Born-Markov 
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approach breaks down. 

4. Discussion and conclusions 

We now turn to consider the imphcations of our results for experiments on 
nanomechanical-SET systems and to explore the analogies which exist with other 
physical systems. 

In an experiment the overall state of the resonator would be determined by the 
combined effects of the SSET and the external surroundings of the resonator other than 
the SSET, which can be modelled by a damping rate, 7e, and a temperature, Tg. The 
overall effective temperature and damping rate of the resonator would be given by the 
weighted averages [S], 

lesTeS = TSET^SET + 7e^e (25) 
7eff = 7SET + 7e- (26) 

From these relations, it is clear that when 7set becomes negative 7efr will decrease, 
eventually leading to a dynamic instability in the state of the resonator when < 0. 
However, we should point out that the model described here is not sufficiently robust 
to explore precisely what would happen in this regime. 

The expressions for the effective temperature and damping rate also make clear that 
for 7sET ^ 7e, the effective temperature of the resonator will be dominated by Tset- 
This has a number of important implications given the minima in TgET which develop 
close to the JQP and D JQP resonances. For example, the minimum value of Tset which 
is simply proportional to the quasiparticle decay rate can easily be as low as ~ lOOmK, 
suggesting that the electronic back-action on the resonator could be demonstrated in a 
dramatic way by using the SSET to cool the resonator when Tg > Tset- Indeed a close 
analogy can be made between the temperature of a resonator coupled to a SSET and 
that of (two-level) atoms undergoing Doppler cooling due to counter-propagating laser 
beams [23 123- 

In Doppler cooling, the interaction between laser beams and two of the atoms' 
energy levels leads to an effective damping of the atomic translational motion. The 
atoms are cooled to a temperature that depends on the detuning of the laser light from 
the atomic resonance||. A, and the decay rate of the excited atomic state, Fg, given by 
the relation 

This equation has an almost identical form to that for Tset close to the JQP resonance 
[equation (j^ ]. It is interesting to note that a more direct analogy between Doppler 
cooling and a system consisting of a resonator coupled to a Cooper pair box (CPB) [21] 
addressed by an additional, fixed voltage gate was suggested recently jHO] • The effective 
temperature of the resonator in that case takes exactly the same form as that near the 

II Note we define the detuning, A, with the opposite sign to that given in Ref. |23| . 
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JQP resonance when an appropriate ac- voltage is applied to the extra fixed gate jHU] , 
although the relevant decay rate is not related to quasiparticle tunnelling and the 
associated expression for the damping of the resonator takes a different form to that 
considered here. 

Another system that might be expected to have very similar dynamics to the 
S SET- resonator system is that of a double quantum dot (DQD) gated by a mechanical 
resonator jST]. Since DQDs display an electronic resonance which is in many ways 
analogous to the JQP cycle, the dynamics of a resonator coupled to a DQD should be 
very similar to that of one coupled to a SSET near the JQP resonance. Experiments on 
DQD-resonator systems have not yet been performed, but work on gating individual 
quantum dots by mechanical resonators [32] demonstrates that such systems are 
certainly a feasible prospect. 

Our work raises a number of interesting questions for future research. For example, 
in the regime that we examined, the strongly non-linear properties of the quasiparticle 
tunnelling rates played no role. It would be interesting to examine the dynamics 
of the resonator for parameters close to the threshold for quasiparticle tunnelling so 
that the motion of the resonator itself could control whether or not tunnelling occurs. 
Also of interest is our finding that the self-consistent Born-Markov approximation 
method apparently cannot capture the physics of the system when it is tuned extremely 
close to the JQP resonance; this raises the intriguing possibility that the resonator 
dynamics in this regime differs substantially from the thermalized dynamics described 
here. Furthermore, we have not yet examined the conditions under which Born-Markov 
approximation breaks down for the DJQP resonance, though this is an area we plan to 
explore in future work. 

In conclusion, we have found that, like a normal state SET, a superconducting SET 
in the vicinity of either the JQP or DJQP resonance acts on a nanomechanical resonator 
like an equilibrium thermal bath for sufficiently weak electro-mechanical coupling and for 
sufficiently large separation between the electrical and mechanical time-scales. However, 
the effective temperature, damping and frequency shift of the resonator due to a SSET 
close to the JQP and DJQP resonances take very different forms, both from each other 
and from the normal state SET. In particular, the magnitude and even the sign of the 
effective temperature and damping for the SSET depend very sensitively on where it is 
tuned with respect to either the JQP or DJQP resonance. 
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Appendix A. Master Equations for the DJQP Resonance 



In this Appendix we present the semiclassical master equations for the SSET-resonator 
system in the vicinity of the DJQP resonance. As shown schematically in figure Eb, 
the DJQP cycle [1^1 1201 (HHI involves resonant Cooper pair tunnelling at each junction 
alternating in turn with two quasiparticle tunnelling events. 

The semiclassical master equations for the SSET-resonator near the DJQP 
resonance are derived using the same procedure as for the JQP resonance. Written 
in dimensionless notation, they take the form 
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The relevant energy differences for the quasiparticle and Cooper pair tunnelling 
processes involved in the DJQP cycle are given by 



and 
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respectively. The DJQP resonance occurs when En^2,n = -E'Ar-i,Af+i = 0, i.e. when 
eVds = 2Ec and A^^ = + 1. Notice that the cycle also requires that the quasiparticle 
tunnelling processes are not energetically forbidden, which implies that Ec > 2 A/3, if 
we assume for simplicity that the superconductors are at zero temperature. 

These master equations give rise to a Fokker-Planck equation with the same form 
as Equation (UH), but with different expressions for the coefficients Cf^, C|^, and C|^. 
These coefficients can be calculated numerically, in the same way as those for the JQP 
resonance, to give the results presented in section 3.3. 
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